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Why is the ground state electron configuration for Lithium ls 2 2s ? 
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Abstract - The electronic ground state for Lithium is ls 2 2s, and not ls 2 2p. The traditional 
argument for why this is so is based on a screening argument that claims that the 2p electron 
is better shielded by the Is electrons, and therefore higher in energy then the configuration that 
includes the 2s electron. We show that this argument is flawed, and in fact the actual reason for 
i-C . the ordering is because the electron-electron interaction energy is higher for the 2p — Is repulsion 

than it is for the 2s — Is repulsion. 
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At the heart of our understanding of atoms, molecules, 
and solids is an interplay between electron-ion, electron- 
electron, and ion-ion interactions. Bound states exist in 
spite of repulsive ion-ion interactions so it is really the first 
two interactions that collaborate to produce conglomerate 
forms of matter. [1] If we confine our discussion to atoms, 
then it is clear from the structure of the periodic table 
that a significant (and appropriate) bias has been applied 
in favour of the electron-ion interactions. [2] Historically 
this is no doubt due to two reasons: first, the single elec- 
tron problem, where electron-electron interactions are ab- 
sent, has an exact solution (so it is easy to 'conceptualize' 
bound states of many-electron atoms in terms of these 
solutions), and second, as the atomic number Z of the nu- 
cleus increases, the energy associated with the electron-ion 
interaction (and with the electron kinetic energy) grows as 
Z 2 , and is therefore dominant. Thus, for the most part, 
the characteristics of electronic properties of the elements 
of the periodic table are described in terms of "Hydrogenic 
states," i.e. states in which electron-electron interactions 
are absent. 

Nonetheless electron-electron interactions can have an 
important effect on the properties of elements, particularly 
when they are introduced into a solid state environment. 
[3] In this paper we wish to re-examine a particularly sim- 
ple case, that of the ground state electron configuration 
in Lithium. The electrons in Lithium are said to reside in 
the Is 2 2s (i.e. 2 S) state, which is standard short- hand for 
the antisymmctrizcd state that is more fully described be- 
low. In the absence of electron-electron interactions, the 
ls 2 2s and ls 2 2p (i.e. 2 P) states for Lithium are degen- 
erate. An accurate computation, utilizing Hylleraas basis 



states, [4] makes it clear that when interactions are fully 
taken into account, the 2 S state indeed has a lower energy, 
as is clearly observed. 

The standard argument for why this is the case pro- 
ceeds as follows. [5] The higher angular momentum of the 
2p electron in the 2 P state tends to "push" the 2p electron 
further out from the nucleus than is the case with the 2s 
electron in the 2 5* state. This is based on a semi-classical 
explanation that invokes the centrifugal force associated 
with the higher angular momentum of the electron. Two 
opposing effects now arise. Being further from the cen- 
tre, the 2p electron interacts less with the inner Is elec- 
trons; this lowers the energy. However, the 2p electron is 
also screened more from the nucleus, and cannot take full 
advantage of the electron-ion interaction; this therefore 
raises the energy. The effects of screening are argued to 
outweigh the lessened interaction with the inner electrons, 
and therefore, the 2 S is the preferred state. 

While both effects are rooted in the electron-electron 
interaction, the screening argument puts emphasis on 
the energy lowered through the electron-ion interaction, 
while the opposing argument clearly emphasizes the direct 
electron-electron repulsion between the various electrons. 
[B] In particular, for a small number of particles, and when 
distances are not large, the distinction between screening 
and electron-electron interactions is not clear. In this pa- 
per, we address this issue and reexamine the traditional 
argument for Lithium's preferred state. We will make the 
claim that it is flawed; contrary to the standard argument 
espoused in the previous paragraph, screening is somewhat 
irrelevant and it is the direct electron-electron repulsion 
between electrons that plays a central role in deciding the 
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ground state. 

The first problem with the standard argument is the 
claim that "the 2p electron is further out from the nucleus 
than the 2s electron." This is incorrect. It is true that the 
2p state has a node at the origin, but Kramer's relation 
[7] allows one to compute the expectation values 
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From either it follows that the expectation value of the 
radial extent of the wave function is larger for the 2s than 
for the 2p state, contrary to the semi-classical explana- 
tion. One can argue that expectation values of the inverse 
powers of the radial coordinate are more important for the 
energetics. Then < 1/r >= l/(n 2 a) doesn't discriminate 
between the two, while < 1/r 2 >= 
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indicate that the 2p state is more extended. 

To settle this argument in an unbiased fashion, we follow 
Ref . [ [5]] , and solve for the radial wave function of a parti- 
cle in an effective potential (Coulomb + centrifugal term) , 
using the one dimensional radial equation. The method of 
solution presented in the preceding reference requires the 
use of an infinite square well potential, extending from 
r = to r = a wa n, to 'embed' the potential of interest. 
For sufficiently large o wa ii the low energy eigenvalues and 
eigenstates are insensitive to the existence of this outer 
wall. [S] However, in the present work, in order to see 
which of the two states has larger extension, as measured 
by the energy, we simply have to lower a wa n and monitor 
which of the two energies, E 2s or E 2p is affected more. 
Fig. 1 shows the result of such a calculation, for both the 
2s and 2p, and for the 3s, 3p, and 3d states. There is no 
question that in either case the energy of the s state is the 
most affected by the presence of the embedding potential, 
and hence the most extended. This is in agreement with 
the expectation based on Eq. ([1]) or ([2]) . In this sense the 
2s state in Lithium is more susceptible than the 2p state 
to being screened by the inner Is electrons. Thus, if the 
naive argument advanced above were correct, the ground 
state would be primarily ls 2 2p, not as observed. We now 
examine the interactions of the 'outer' n = 2 electron with 
the Is electrons. 

We first use a simple variational wave function to calcu- 
late the energies of the ls 2 2s and \s 2 2p states in Lithium. 
Our trial wavefunction has the usual determinant form ap- 
propriate to the case where we ignore interactions. It is 
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Fig. 1: (Colour online) Energy of the 2s vs. the 2p states 
(lower two curves) and 3s vs. 3p vs. 3d states (upper three 
curves) calculated as a function of the width of an embedding 
infinite square well potential. For large width, the states are 
degenerate as expected for the Coulomb potential, but as the 
width is lowered, the first state affected in either case is the s- 
state, indicating that it is the more extended of the degenerate 
set. 
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Here a is the Bohr radius, and a(r) and /3(r) represent 
the usual spin states. For Lithium, z\ — z 2 — 3; here they 
are retained as variational parameters. We use this trial 
wavefunction to place an upper limit on the energy of each 
state. The nonrclativistic Hamiltonian for the Lithium 
atom is 
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where the terms corresponding to each part of the Hamil- 
tonian should be clear. Note that Z is the nuclear charge. 
Contributions to the energy are divided into the follow- 
ing categories: K =< i? k inotic >, ^nuclear =< -^nuclear >, 
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and Uuu, U ls2s , U ls2p , Ex ls2s , and Ex ls2p refer to the di- 
rect Coulomb and exchange energies for the orbitals listed. 
Then the total energies for each state are 



Es — K + U nuc \ eai + Uisis 

Ep = K + ^nuclear + Crisis 



2J7i s2s + Ex\ s2s (7) 
2U ls2p + Ex ls2p . (8) 

Calculations are straightforward; with Eq = h 2 /(2ma 2 ) 
(= 1 Rydberg), the three terms in common for the two 
states are 
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Then for the two trial states, we have 
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respectively. The exchange terms (Ex\ s2s and Ex\ s2p ) 
arise from consideration of the full antisymmetric wave- 
function. For the two trial states, the exchange terms are 
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respectively. In addition, the first state is not normalized. 
Therefore 
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Term 



Perturbative (eV) 
ls 2 2s ls 2 2p 



Variational (eV) 
ls 2 2s ls 2 2p 



K 


275.5 


275.5 


212.1 


200.0 


v nuclear 


-551.0 


-551.0 


-486.7 


-459.9 


Ulsls 


51.0 


51.0 


46.6 


45.7 


U\ s 2s 


17.1 




11.4 




U\ s2p 




19.8 




7.1 


Exch ls2s 


-1.8 




3.9 




Exch ls2p 




-1.4 




-0.1 


<*l*> 


1 


1 


0.98 


1 


Total 


-192.0 


-186.3 


-201.2 


-200.0 



Table 1: Energies found for each term using a perturbative 
calculation (z\ — zi = 3) and a variational calculation using 
Eqs. © and ([SJ. In both cases Z = 3. The optimal z\ and 
Z2 values from the variational calculation are given in the text. 
The actual values for the ls 2 2s and ls 2 2p energies are —203.5 
eV and —201.6 eV respectively. [4] 



The trial wave function remains normalized for the 2 P 
state, i.e. (*p|*p) = 1. 

If the variational parameters are fixed at z\ = z 2 = 
3, the calculation may be considered as a perturbative 
calculation, where the electron-electron interactions act 
as the pertubation; these energies are listed in Tabic [TJ 

More accurate energy estimates may be found if Z\ and 
z 2 are taken as variational parameters. Minimization of 
Eq. yields z x « 2.68 and z 2 « 1.87, while for the 2 P 
state given by Eq. ([5]) we find z\ w 2.69 and z 2 w 1.05. 
The proximity of the last number to Z — 2 = 1, and the 
relative accuracy of the variational calculation (differences 
of 1.1% and 0.8% for the two states) no doubt combine to 
inspire the traditional screening argument. 

However, these variational calculations implicitly ac- 
count for electron configurations beyond the n = 1 and 
n = 2 states under consideration here, More straightfor- 
wardly, inspection of the perturbation-like calculations in 
Table 1 indicates that the electron-electron interaction en- 
ergy between the Is and 2p states (2(Ui s2p + Exi s2p ) w 
36.8 eV) is considerably higher than that between the Is 
and 2s states {2{U\ s2s + Ex\ s2s ) w 30.6 eV). When one 
goes beyond perturbation theory, both wave functions ex- 
pand by incorporating higher energy configurations; these 
serve the purpose of allowing the electrons to avoid one 
another more efficiently, and in so doing, the Is electrons 
effectively 'screen' the 2s or 2p electron. So in an attempt 
to ameliorate the (rather large) electron-electron energy 
cost for the 2p electron, the variational calculation simply 
places it further from the nucleus, and succeeds in low- 
ering this energy, but at a cost of what we have called 
^nuclear (—551.0 gets raised to —459.9). The same thing 
happens for the 2s electron, but to a lesser extent, in part 
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because the 2s electron is already further extended than 
the 2p electron. Nonetheless, this does not succeed in re- 
versing what was evident from the first two columns - 
the ls 2 2p state simply has a higher electron-electron en- 
ergy compared to the ls 2 2s state, and therefore the latter 
has lower energy. 

Phrased this way, it should be evident that 'screening' 
is the 'sacrifice' one makes to lower the overall energy. In 
both the 2s and the 2p case, lowering the overall energy 
was accompanied by screening, i.e. a raising of {/nuclear 
(while the overall energy decreased). This is qualitatively 
different than the standard argument that says that the 2p 
state is more screened to begin with (false!) and therefore 
its energy is higher. 

In summary, we have argued that the ground state elec- 
tron configuration for Lithium is ls 2 2s, and not \s 2 2p in 
spite of screening, not because of screening. Expansion of 
an electron cloud surrounding a nucleus generally occurs in 
an effort to minimize the Coulomb repulsion experienced 
by these electrons. [TU] It follows that the more repulsive 
these interactions are, the more will the electrons expand 
away from one another. In so doing, they also become 
more efficiently screened by the remaining electrons; this 
is a sacrifice that is made in spite of the electron-ion en- 
ergy increase, because the total is thereby decreased. The 
original and incorrect argument stated that the Hydro- 
genic 2p state is more screened than the Hydrogcnic 2s 
state, and this is why the latter results in a lower energy. 
Our argument clearly attributes the deciding factor to the 
electron-electron repulsion, and increased screening is a 
necessary (energy raising) by-product that is nonetheless 
tolerated by the ground state. 

* * * 

This work was supported in part by the Natural Sciences 
and Engineering Research Council of Canada (NSERC), 
and by the Teaching and Learning Enhancement Fund 
(TLEF) at the University of Alberta. 

REFERENCES 

[1] We single out the ion(s) in this discussion because they are 
much more massive than the electrons; in other words we 
assume a Born-Oppenheimer approximation. 

[2] Note that in the case of solids the dilemma of focussing on 
electron-ion vs. electron-electron interactions faced Bloch 
in his initial formulation of "Bloch's Theorem;" he focussed 
on the former basically because it was "the simplest as- 
sumption." See F. Bloch, Proc. R. Soc. Lond. A371, 24 
(1980), and, for a more recent discussion, J.E. Hirsch, Int. 
J. Mod. Phys. B23, 3035 (2009). 

[3] This statement is meant to have more content than might 
appear at first glance; it goes without saying that electron- 
electron interactions result in subtle correlations in solids 
from which new phases of matter, such as magnetism or 
superconductivity may result. What we mean here, how- 
ever, is that the atomic (i.e. local) properties of electronic 



states can be significantly altered by electron-electron in- 
teractions. An example is the necessity of acounting for 
expanded atomic orbitals, when those orbitals are doubly 
occupied, as first suggested clearly by J.E. Hirsch, Phys. 
Rev. Lett. 87, 206402 (2001). 

[4] See, for example, Z.-C. Yan and G.W.F. Drake, Phys. 
Rev. A52, 3711 (1995), and more recently L.M. Wang, Z.- 
C. Yan, H.X. Qiao, and G.W.F. Drake, Phys. Rev. A85, 
052513 (2012), and references therein. 

[5] This argument is given in many standard texts, for exam- 
ple, D. J. Griffiths, Introduction to Quantum Mechanics 
(Pearson/Prentice Hall, Upper Saddle River, NJ, 2005), 
2nd edition (see page 215). 

[6] Note that in the case of Helium, while the second elec- 
tron is primarily in the Is state because that minimizes 
the electron-ion energy, there is nonetheless almost 10% 
probability that either electron is not in the Is level, specif- 
ically because this lowers electron-electron interaction en- 
ergy. Thus in the case of Helium, screening (i.e. some of 
the electron cloud is pushed further out) occurs to lower 
the energy, not the other way around. 

[7] See Ref. [ [5]], p. 289, problems 6.34 and 6.35, where these 
relations are derived. 

[8] B.A. Jugdutt and F. Marsiglio, "Solving for 3D potentials 
using matrix mechanics with a truncated basis set," un- 
published. 

[9] The existence of the wall at r = does not effect the solu- 
tion since the radial wave function u(r) = rR(r) is required 
to be zero there anyways. 

[10] See the reference in Ref. [ [3]], along with J.E. Hirsch, 
Phys. Rev. B65, 184502 (2002). 



p-4 



